Time reversal of pseudo-spin 1/2 degrees of freedom 
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^Abstract 

^ We show that pseudo-spin 1/2 degrees of freedom can be categorized in two types according to their behavior under time reversal. 
One type exhibits the properties of ordinary spin whose three Cartesian components are all odd under time reversal. For the second 
\^ type, only one of the components is odd while the other two are even. We discuss several physical examples for this second type of 
PsJ pseudo-spin and highlight observable consequences that can be used to distinguish it from ordinary spin. 
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1. Introduction 

The behavior under time reversal (TR; also denoted rever- 
sal of motion) is one of the most fundamental characteristics 
of a quantum system LI, 2J. It determines, e.g., level degen- 
eracies ^ and the statistics of energy-level spacings in 
closed systems, electric transport in phase-coherent quantum 
circuits and the possible channels for pairing of electrons 

to form a superconducting Cooper-pair condensate |7]. For- 
mally, the TR operation can be represented by an anti-unitary 
operator that is, however, specific to the particular choice of 
base kets |2, 3]. Quite generally, we may write 



e^uc 



(1) 



where U denotes a unitary operator and C is complex conjuga- 
tion. For simple quantum systems, the explicit form of the TR 
operator 9 is well known [2, Recent efforts were aimed at 
generalizing the TR operation to more complex systems, e.g., 
those having internal degrees of freedom |8]. 

We focus here on TR of quantum systems that carry an ef- 
fective SU(2) (pseudo or real) spin degree of freedom that may 
be half-integer or integer. We show that SU(2) symmetry al- 
lows for two fundamentally different behaviors under TR. As 
TR is an antiunitary symmetry independent of the unitary sym- 
metry elements in SU(2), it needs to be determined based on 
physical considerations which TR behavior applies to a partic- 
ular system. This result has important consequences for effec- 
tive pseudo-spin descriptions that are widely utilized. Classic 
examples include Schwinger's oscillator model of angular mo- 
mentum 191] (see also Sec. 3.8 in Ref. lH), nuclear isospin |10], 
and the ammonia molecule [IIU. More recently, the pseudo- 
spin concept has been ubiquitous in the context of quantum 



information processing (l^. Other pseudo-spin-carrying enti- 
ties of current interest include the massless Dirac-electron-like 
quasiparticles in graphene ifisll and the persistent spin helix in 
quasi-twodimensional semiconductor systems with fine-tuned 
spin-orbit couplings fU, Ts']. We will discuss these particu- 
lar examples and elucidate experimentally observable ramifica- 
tions for the two different types of pseudo-spins. 

2. Time reversal of a pseudo-spin: General properties 

We start by considering the textbook example of an SU(2) 
angular-momentum algebra involving the three operators Jj 
with j = 1, 2, 3, satisfying the commutation relations 



[jj, Jii] - iejkiJi 



(2) 
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In general, kinematically relevant physical quantities are either 
even or odd under TR 1,16.1 . Allowing for both possibilities 
for each operator Jj, we write 6Jj0'^ - ^jJj, with - +1. 
As the commutators (|2]l need to be preserved under TR, the 
three coefficients cannot be independent. Rather, they must 
satisfy the condition ^1^2 - -^3, which restricts the possi- 
ble TR behavior to two cases: (f) - ^2 - ^3 = -1, or 
(b) = -^2 = -^3 = -1, with permutations of indices al- 
lowed. Case (f) impUes that all operators Jj are odd under TR, 
which is the behavior found for the Cartesian components of 
orbital and ordinary-spin angular momentum |3]. In case (b), 
only one of the operators Jj is odd under TR (without loss of 
generality chosen here to be Ji), and the other two {J2 and 73) 
are both even. For an SU(2) invariant system of type (b), one of 
the Cartesian components of J is thus always distinguished by 
its behavior under TR. 

The TR behavior associated with cases (f) and (b) generally 
leads to qualitatively different physical properties. For half- 
integer (pseudo-) spin systems, a basic feature distinguishing 
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the two cases is given by the fact that 0^ - - I, whereas case (b) 
implies 0^ - +1 UtIi . To prove these relations, we construct 
the TR operator for cases (f) and (b) in the usual representa- 
tions jsl], where the matrix elements of /] and /} are real, while 
those for J2 are imaginary. Given the general relation expressed 
in Eq. ([TJ, it is only necessary to find unitary transformations 
Uf and Ub (i.e., rotations in spin space) that yield the required 
transformation properties of the operators Jj. As J2 is odd un- 
der complex conjugation while Ji and J3 are even, it is straight- 
forward to find Uf = exp(/;7r72) and Uh - exp{inJ3). Noting 
that COfC- Uf, we find 0^ = 0^= exp {linJ^) amounts to a 
2n rotation in spin space, which yields 0^ - -\ for half-integer 
spin. Conversely, CtJhC -U^^ and thus 0^ - tJhU^^ = 

The sign of r has direct experimental consequences for 
the quantum interference of amplitudes for TR-related physical 
scenarios. Such interference plays a role, e.g., in particle inter- 
ferometers 1 1 8, 19 1 and for the phase-coherent propagation of 
electrons and light through disordered media 1^. 



3. Physical realizations of case-(b) pseudo-spin 

Whether a (pseudo-) spin degree of freedom belongs to case 
(f) or (b) depends on the details of the particular system under 
consideration. It is well-known |3] that orbital and ordinary- 
spin angular momentum belong to case (f). In contrast, the 
possibility to have the TR properties of case (b) has been under- 
appreciated, even though physical realizations exist. For exam- 
ple, it is well-known that the third Cartesian component of nu- 
clear isospin must be even under TR by virtue of electric-charge 
conservation, and some implications have been considered in 
previous work 10,01 • Below we further elucidate the general TR 
properties of type-(b) pseudo-spins using isospin as an example, 
and also Schwinger's bosonic model of angular momentum dgj, 
which is equivalent to the isotropic two-dimensional (2D) har- 
monic oscillator 1 2^ 21 1 . Finally, we show that case (b) applies 
also to the persistent spin hehx lll411 . 



3.1. Nuclear isobaric spin 

The isospin- 1/2 model is used to describe states of the nu- 
cleon LLQJ, with eigenstates of ^3 being associated with the pro- 
ton and neutron, respectively. The electromagnetic interaction 
distinguishes between the two states, and conservation of elec- 
tric charge requires that all relevant states are eigenstates of ^3. 
In addition, charge conservation mandates the invariance of J^, 
under TR. In the usual representation where is real, this im- 
plies J-i 0"' - U J3 U^^ - 73, i.e., U must be a rotation around 
the 3-axis in isospin space. As the physically relevant isospin 
states are eigenstates of 73, U was not specified further in pre- 
vious work [2]. However, it is straightforward to verify that 

^ uu-^ ^ +1, i.e., isospin is an example for the type-(b) 
pseudo-spin. 

3.2. Isotropic 2D harmonic oscillator and Schwinger model 

It was realized early on ll20ll that the three dynamic invari- 
ants of the isotropic 2D harmonic oscillator 12111 correspond 



to a dynamical SU(2) symmetry. Using dimensionless coordi- 
nate and momentum variables, in which the Hamiltonian reads 
H = (pj + q\)l2 + (P2 + q\)l'2, the three conserved quantities 
can be expressed as /i = {.qxPi-qiPiM^-, J2 - {qiq2 + PiPi)/'^, 
and 73 = [p2 + q\ — {p\ + q\)\ /4. They correspond to the orbital 
angular momentum of the oscillator (7i), the energy diff'erence 
for motions in the two perpendicular in-plane directions (73), 
and the phase difference between oscillations in those directions 
(72). Straightforward calculation based on the canonical com- 
mutation relations {qj,pk\ - iSjk and {qi,qj\ - {pi,Pj\ - 
establishes that the Jj satisfy Eq. Furthermore, the equiv- 
alence between this system and Schwinger's bosonic model of 
angular momentum |3, 9] becomes apparent when the quanti- 
ties J j are expressed in terms of creation and annihilation oper- 
ators fl^ - (qk - ipk)l V2, fl* = {qk + ipk)l V2 for the two ID 
oscillators. 

As the coordinate (momentum) components are even (odd) 
under TR, it follows that J\ is odd but both 72 and 73 are even 
so that this system is a realization of type (b). Even when the 
Schwinger-model description is applied to a general two-level 
system and, thus, the underlying oscillator degree of freedom 
is abstract, the definition of 7? = (a2'3'2 - a\a\)l2 in terms of 
the occupation numbers of the two levels implies that 73 must 
be even under TR. Hence, many pseudo-spin models used in 
condensed-matter physics ll22ll and quantum optics ll23ll belong 
to type (b). 

3.3. Persistent spin helix 

The persistent spin helix is a recently discovered fl4, f51 
collective excitation present in certain quasi-2D semiconduc- 
tor systems with fine-tuned spin-orbit couplings. In terms of 
second-quantized operators ^^.j^^ and cii.kj,) that, respectively, 
create and annihilate an electron with wave vector k and spin-up 
(spin-down), the following operators associated with the persis- 
tent spin helix are defined 

= I]^k+Q,T^k,i, 5q = J]c;^Ck+Q,T , (3a) 

k k 



7 ^VkA^^'^'^Ki^^-^ ' 
^ k 



(3b) 



where is a function of the spin-orbit coupling strength in the 
system 1141. If we define S ^ - (Sq + Sq)/2 and 5, - (5i- 
5q)/(2/), it can be shown that the components S j obey Eq. ©. 
Using the fact that both k and spin t, I are odd under TR, we 
find immediately §5; 0"' - -S~. Similarly, we get 05 q = 
Sq which implies that S x and 5^ are even under TR. Hence, 
the SU(2) degree of freedom associated with the persistent spin 
helix is a type-(b) pseudo-spin. 

4. Conventional and unconventional TR of 2D massless 
Dirac particles 

In the previous section we discussed examples, where the 
components Jj of the pseudo-spin were conserved, [H, Jj] - 
so that the TR properties of Jj could be discussed separately 
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Figure 2: (a) Honeycomb lattice of graphene. Atoms in sublattice A (B) are 
marked witli open (closed) circles, (b) Brillouin zone and its two inequivalent 
comer points K and K'. The remaining comers are related with K or K' by 
reciprocal lattice vectors, (c) Dispersion E(k) near the K point. 



Figure I: (a) Bloch sphere for the pseudospin of 2D massless Dirac particles. 
The pseudo-spin for wave vector k is characterized by a rotated coordinate 
system ((Tk, <t±, <t-) with cTk parallel to k. (b) Time reversal (TR) of pseudospin 
parallel to k can be related to TR of pseudospin parallel to the x axis via two 
rotations yi and 73 around the z axis by angles -^Ok and </>k. Depending on 
whether TR of pseudospin parallel to the x axis is achieved via a rr rotation 
around the y axis (7') or around the z axis (y^). one obtains &f or dj,. 



from the orbital motion. As a classic example where the 
(pseudo-) spin degree of freedom is coupled to the orbital mo- 
tion, we now discuss particles confined to the xy plane that are 
described by a massless Dirac Hamiltonian. Using a plane- 
wave (k) basis, the Hamiltonian becomes H{k) = tivkcTk, with 
velocity v and the (pseudo-) spin operator ixk 



(Tk = (Ti- COS ip\i + cr,. sm ^k = 











(4) 



Here the ctj denote the familiar Pauli matrices, and i^k is the 
angle between k and the k^- axis [see Fig.lita)]. The eigenvalues 
of H{k) are +hvk, and the corresponding eigenstates are 



|k,±) = e 



ikri 



V2 



>-'Vk/2 



±e' 



'iPk/2 



(5) 



It is reasonable to enable a unique identification of the eigen- 
states by restricting -n < (f]^ < n. 

Reversal of motion should map the eigenstates of H(k) as 
follows: 

^|k,±) = /7(k)|-k,±) , (6) 

where //(k) stands for an arbitrary phase, i.e., h should reverse 
the wave vector while preserving the energy. Inspection of 
Eq. (|5]l yields the relation |±)-k = ' sgn(^k) l+)k, consistent with 
cr_k = -CTk. Thus (Tk needs to be odd under the TR operation, 
to ensure TR invariance of i/(k). For a particle confined to two 
spatial dimensions, this condition does not uniquely specify the 
TR operation. Rather, two possible scenarios exist for how the 
direction of |k, ±) can be reversed for any given k. This is illus- 
trated in Fig. [lib)- Mathematically, the two different TR opera- 
tions are given by 

= exp(- 5(^k o-J exp( \-k cr,,) C exp( ^(^k cr^) 

= exp(57rcr,,)C = /cr^^ C , (7a) 

= exp(-^(^k o-;) exp[^ sgn(v3k)7ro-J C exp(^^k cr-^ 

= exp(-^[^_k + i^k] o"z) C . (7b) 

It is straightforward to verify that iy(k) is invariant under 
both types of TR transformation. However, )?f and ^b can be 



distinguished by the way the pseudo-spin operators cr- and cr^ 
are transformed, [ctj^ is the pseudo-spin in-plane component 
perpendicular to o"k, see Fig. [Ha).] More specifically, both ctj. 
and cTj^ turn out to be odd under )?t, whereas they are even 
under ^b- Hence, which of the two operators )^b,f constitutes 
the proper TR operation for a particular system will depend on 
whether the actual physical quantities represented by cr, and ctj^ 
are even or odd under TR. 

Before closing this Section, we briefly discuss time-reversal 
properties of the sublattice-related pseudo-spin degree of free- 
dom carried by electronic quasiparticles in graphene 1 13]. This 
material consists of a sheet of carbon atoms arranged on a 
honeycomb lattice, with two inequivalent sublattices labeled A 
and B [Fig. 12 a)]. Its Brillouin zone is also hexagonal, with two 
inequivalent corner points K and K' - -K [Fig.|2jb)]. Most 
remarkably, the energy dispersion near K and K' is linear, and 
the Fermi energy of undoped graphene is exactly at the crossing 
point [Fig. |2tc)]. Hence, the low-energy electronic excitations 
in graphene are considered to be analogs of 2D massless Dirac 
electrons ElMSj. 

In situations where Umklapp processes are absent, it is 
tempting to treat the 2D Dirac cones at the K and K' valleys 
separately and also define an effective intra-valley time-reversal 
operation l27l 12811 . As the pseudo-spin operator cr, represents 
the quasiparticle location on sublattice A or B, it must be even 
under time reversal (as expected for a real-space position oper- 
ator). i?b of Eq. ( iTbl i exhibits the required property to leave cr^ 
invariant (in contrast to i?f ), suggesting it to be a proper intra- 
valley time-reversal operator However, it should be noted that 
the ordinary TR for Bloch electrons in graphene is given by 
- C which couples the two valleys Il28i1 . In order to apply 
the discussion of this Section to graphene, a restriction of anti- 
unitary operators to individual valleys needs to be formulated. 

5. Observable signatures of the two types of pseudo-spin 

Whether a particular realisation of pseudo-spin belongs to 
type (f) or (b) is not a purely academic question, as their differ- 
ent TR properties affect physical observables. For example, in 
the case of 2D Dirac particles considered in the previous Sec- 
tion, a (possibly random) potential y(r)crj. does not break ^b- 
TR symmetry, whereas such a potential will break the df-TR 
symmetry f2^. As a result, the level statistics of chaotic quan- 
tum dots formed by a non-integrable mass-confinement of 2D 
Dirac particles will be different in the two cases ll29ll . While 
this example is specific to the case of 2D Dirac particles, we 
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will now discuss a more general observable difference exhibited 
by half-integer pseudo-spin particles of type-(f ) and (b) that are 
scattered by a spin-independent random potential but are free 
otherwise. The quantum interference of TR-related backscatter- 
ing amplitudes turns out to be different for the two pseudo-spin 
types and can thus serve to distinguish them experimentally. 

We denote scattering states of the particle by |k, s), where 
k denotes the wave vector and s is the pseudo-spin quantum 
number w.r.t. some specific basis. The probability amplitude for 
a particular backscattering process involving n scattering events 
due to a disorder potential is then proportional to 

■^ki.si,...,K,s„ = <-k, i|k„, i„)<k„,i„|k„_i,i„_i)... 

...<k2,52|ki,5iXki,5i|k,i). (8) 

In the typical situation where disorder potentials are invariant 
under TR, each such process has a "partner" process 

^k„..,,...,k„..„ = <-k,i|0(ki,ii)X0(ki,ii)|§(k2,i2)>... 

. . . <0(k„_i, i„_i)|§(k„, s„)}{6{k„, s„)\k, s) , (9) 

where scattering occurs in time-reversed order. Using the rela- 
tions H {a\b} = {6b\6a} and |0(k, s)) = |-k, s), we find 

^k,,.„...,k„A, = <-k, 5|02(k„, i„)X§2(k„, 5„)| . . . 

...|§2(k,,i,)X0'(kl,5i#(k, 5)> 

= sgn(02)y[k,_,_ . (10) 

The total probability for backscattering of a particle is the 
modulus square of the sum over the probability amplitudes of 
all possible back-scattering processes. If the sign of 0^ is posi- 
tive (negative), scattering processes related by TR will interfere 
constructively (destructively), leading to reduced (enhanced) 
particle transmission through the medium. Thus we suggest 
that quantum-coherent transport provides an avenue for dis- 
tinguishing between half-integer type-(f) and type-(b) pseudo- 
spins, because 6^ = -0^ = -1. Direct observation of the 
backscattering probability would facilitate such an experimen- 
tal distinction. 

6. Conclusions 

We have shown that (pseudo-) spin degrees of freedom can be 
classed into two types according to their time-reversal proper- 
ties. For type (f), which includes ordinary spin angular momen- 
tum, all three spin components are odd under time reversal. For 
type (b), only one of the components is odd under time reversal 
while the remaining two are even. This case includes nuclear 
isospin and the dynamic SU(2) symmetry of the 2D isotropic 
harmonic oscillator 

As an example of a system for which (pseudo-) spin and or- 
bital motion are coupled, we discussed time reversal of mass- 
less Dirac particles confined to move in two spatial dimensions. 
It was found that the dynamics of the system again allows for 
time-reversal operations of the two types (f) and (b). Only 
by considering the properties of additional observables repre- 
sented, e.g., by the z component of pseudo-spin is it possible to 
uniquely determine the form of the time-reversal operation. 



Which type of time-reversal operation is realised in a particu- 
lar system has measurable consequences. For example, whether 
itinerant particles with half-integer (pseudo-) spin belong to 
type (f) or (b) manifests itself in an experimentally observ- 
able way by the relative sign of probability amplitudes associ- 
ated with time-reversal-related back-scattering scenarios. This 
and several more observable ramifications associated with ran- 
domness/chaotic dynamics js]] arise because the squared anti- 
unitary operators associated with the time-reversal transforma- 
tion in the two cases differ by their sign. 
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